A renormalizat ion-group analysis of the interacting resonant level model at finite bias: 
Generic analytic study of static properties and quench dynamics 
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Using a real-time renormalization group method we study the minimal model of a quantum dot 
dominated by charge fluctuations, the two-lead interacting resonant level model, at finite bias volt- 
age. We develop a set of RG equations to treat the case of weak and strong charge fluctuations, 
together with the determination of power-law exponents up to second order in the Coulomb in- 
teraction. We derive analytic expressions for the charge susceptibility, the steady-state current 
and the conductance in the situation of arbitrary system parameters, in particular away from the 
particle-hole symmetric point and for asymmetric Coulomb interactions. In the generic asymmetric 
situation we find that power laws can be observed for the current only as function of the level posi- 
tion (gate voltage) but not as function of the voltage. Furthermore, we study the quench dynamics 
after sudden switch-on of the level-lead couplings. The time evolution of the dot occupation and 
current is governed by exponential relaxation accompanied by voltage-dependent oscillations and 
characteristic algebraic decay. 

PACS numbers: 05.60. Gg, 71.10.-w, 73.63.Kv, 76.20.+q 



I. INTRODUCTION 

The standard set-up for a quantum dot consists of a 
small quantum system described by a finite-dimensional 
Hilbcrt space which is coupled to several infinitely large 
reservoirs via energy and/or particle exchange. A dif- 
ference in the chemical potentials of the reservoirs will 
generically lead to particle transport and thus a finite 
current through the dot. Here we will study the arguably 
simplest but non-trivial quantum dot system, namely the 
interacting resonant level model (IRLM) . It is given by a 
local level coupled to two leads of non-interacting spin- 
less fermions. The fermions can hop on and off the level. 
In addition, there is a Coulomb interaction between the 
level and the reservoirs (see Fig. [T}. The IRLM consti- 
tutes the minimal model for a quantum dot dominated 
by charge fluctuations as spin degrees of freedom are not 
taken into account. 

Originally the (one-lead) IRLM was introduced in- 
dependently by Wiegmann and Finkelstein 1 as well as 
Schlottmann 2 to study the anisotropic Kondo model. 
They generalized earlier works by Anderson et al£ at 
the Toulouse points^, where the Coulomb interaction be- 
tween the level and the reservoir vanishes. In particular, 
m Refs. Q and it was shown that the IRLM and the 
anisotropic Kondo model possess the same partition func- 
tion in the so-called long-time approximation and thus 
share the same universal low-temperature characteristics. 
Equilibrium properties like the static and dynamic sus- 
ceptibilities and the relaxation rate of the IRLM have 
been intensively studied in the early 1980's using the 
Bethe Ansatasi as well as renormalization group (RG) 
techniques 6 . The equivalence between the IRLM and the 



anisotropic Kondo model can be shown by bosonization 
and refermionization of the latter model^. 

Recently the interest in the IRLM has been revived 
as a minimal model to describe non-equilibrium trans- 
port through quantum dots. Initialized by the work of 
Mehta and Andrei^ the model has been investigated us- 
ing the Hershfield F-operator— , Keldysh perturbation 10 
and scattering theor y 11 ' 12 , field theory approaches^ 3 - - — , 
the numerical renormalization group method (NRG) 16 
and the time-dependent density matrix renormalization 
group technique (TD-DMRG)^^. Most of these studies 
were performed at the special point of particle-hole and 
left-right symmetry. The quantity of main interest has 
been the steady-state current through the resonant level. 
The main conclusions were that (i) at sufficiently large 
bias voltages a negative differential conductance appears, 
and that (ii) in the scaling limit, where all bare energy 
scales are much smaller than the bandwidth of the leads, 
the current decreases as a power law in the applied volt- 
age. However, only at the self-dual points it has been 
possible to derive closed analytic expressions for the cur- 
rent as a function of the applied voltage. 

Recently, perturbative RG techniques in non- 
equilibriumi^— have been applied to obtain more in- 
sight into the physics of the IRLM at finite bias. In 
Ref. [l6| a poor man scaling analysis has been performed 
up to next-to-leading order providing power-law expo- 
nents up to second order in the Coulomb interaction. 
In a non-equilibrium situation, the RG flow was cut off 
heuristically by the voltage which induced an emergent 
power-law behavior of the current as a function of the 
voltage. Subsequently, this analysis has been put on a 
more firm basis by the application of recently developed 
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RG methods in non-equilibrium, the functional RG— and 
the real-time RG method^, see Ref. [2l| for a recent re- 
view. A short summary of the main results of the two 
methods for the IRLM has been presented in Ref. [ill 
where a leading order expansion has been performed giv- 
ing rise to power-law exponents linear in the Coulomb 
interaction. The results were compared to numerically 
exact NRG and DMRG methods and a good agreement 
has been observed for moderate Coulomb interactions. 
In particular, the conclusion was drawn that power-law 
behavior does not take place in the generic case of asym- 
metric Coulomb interactions between the dot and the 
left and right reservoir. In addition, the scaling behav- 
ior at resonances away from the particle-hole symmetric 
point has been reexamined 9 . Details of the functional 
RG method have been presented in Ref. |rjj 

In this paper we will present an extended version of 
Ref. HH concerning the real-time RG method, supple- 
mented by a generic treatment of strong charge fluctu- 
ations, a next-to-leading order analysis in the Coulomb 
interaction and a new result concerning power laws as 
function of the level position away from the particle-hole 
symmetric point. The real-time renormalization group 
in frequency space (RTRG-FS) 2 — has recently been in- 
troduced in the theory of dissipative quantum systems. 
It provides a powerful tool in the description of non- 
equilibrium transport, in particular the relaxation and 
decoherence rates naturally arise within the proposed for- 
malism. Previous applications to the Kondo mode l 22 i 23 
in the weak coupling regime are here generalized to in- 
clude charge fluctuations in strong coupling. To this end 
we develop RG equations, where we expand all quantities 
around zero Matsubara frequency, in contrast to previous 
treatments 22 ^, where a systematic expansion around 
the poor man scaling solution has been performed. The 
RG equations are set up in a generic form, which can 
also be used for other models with strong charge fluctu- 
ations. In particular, for the IRLM we demonstrate that 
this scheme allows the study of observables close to res- 
onances where the tunneling rate is the only relevant en- 
ergy scale quantifying charge fluctuations. Furthermore, 
we extend the analysis in Ref. by including subleading 
terms, which gives the exponents of power laws consis- 
tently up to second order in the Coulomb interaction. 
A corresponding comparison of the power-law exponent 
with NRG results for the charge susceptibility in equi- 
librium at the particle-hole symmetric point yields excel- 
lent agreement. We present approximate analytical solu- 
tions which are confirmed by numerically integrating the 
corresponding full RG equations and which describe the 
steady state as well as the quench dynamics for arbitrary 
system parameters. Thereby various microscopic cutoff 
scales of the RG flow can be identified, which is essen- 
tial for the precise determination of the scaling behavior 
of observables. In particular, we derive closed analytic 
expressions for the charge susceptibility, the steady-state 
current and the differential conductance. We find (i) a 
negative differential conductance for arbitrary system pa- 
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FIG. 1. (Color online) Sketch of the interacting resonant level 
model. A local level with energy eo is coupled via hoppings 
^l/r an d Coulomb couplings U^j R to two spinless fermionic 
reservoirs held at chemical potentials Hl/r = ±V/2. 



rameters, (ii) that for asymmetric Coulomb interactions 
the current does in general not follow a power law as a 
function of the bias and is recovered only in the limit of 
extremely large voltages, (iii) that at resonance, i.e. when 
the level position is aligned with one of the Fermi levels 
in the leads, the current does not follow a power law even 
in the symmetric model, and (iv) that the current or the 
linear conductance reveals a power law as function of the 
level position in the generic case, i.e. even for asymmet- 
ric Coulomb interactions and/or asymmetric tunne ling 
couplings. The latter result was not reported in Ref. [l8|. 

In addition, we use the analytical solution of the 
RG equations to study the quench dynamics in the 
IRLM, where we assume the couplings to the leads to 
be switched on suddenly. We derive closed integral rep- 
resentations for the resulting time evolution of the dot 
occupation and the current. The most notable charac- 
teristics of the time evolution of both observables are: 
(i) the relaxation towards the stationary values is gov- 
erned by two different decay rates describing the charge 
relaxation on the level and its broadening induced by 
the coupling to the leads, respectively, (ii) the voltage 
appears as an important energy scale for the dynamics 
setting the frequency of an oscillatory behavior, and (iii) 
the exponential decay is accompanied by an algebraic be- 
havior with an interaction-dependent exponent. Similar 
results have been obtained recently for the dynamics of 
the non-equilibrium Kondo model 2 ^, showing that these 
features are generic. 

The paper is organized as follows: In Section |TT] we 
introduce the IRLM and discuss its description in Liou- 
ville space. In Section IIIII we summarize and solve the 
RG equations. In Section IIVI we present the results for 
steady-state quantities as well as for the time evolution. 
Here we also provide a simple derivation of the appear- 
ance of the negative differential conductance. Techni- 
cal details together with the generic derivation of non- 
equilibrium RG equations in the regime of strong charge 
fluctuations are reported in the Appendix. 
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II. MODEL 



The Hamiltonian of the interacting resonant level 
model (IRLM) depicted in Fig. Q] is given by 



H — H res + Hi 



V 



where 



ka 



(eft + p a )a\ a a kc 



(1) 



(2) 



describes a set of semi-infinite fermionic reservoirs with 
chemical potentials p a . In the case of two reservoirs 
a = L/R, we choose Pl/r. = ±V/2. Standard second- 
quantized notation is used, and the energies t ka are re- 
stricted to a finite band of width D. The dot Hamiltonian 
reads 



H D = eo^c 



(3) 



and the fermionic level is coupled to the reservoirs via 



(0) 



> n { °> k 
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tc ~o )Y,^kY, :a L a k> a -- , (4) 
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a Pa kk 

where : . . . : denotes normal-ordering, and and Ua ^ 
are real. In the following we denote the bare parameters 
by the super-script ^ . In the scaling limit the details of 
the frequency dependence of the local density of states 
in the reservoirs p a (oj) do not play a significant role as 
long as it is sufficiently regular on the energy scale of the 
applied voltage, which allows us to appoximate it by a 
constant. Following Ref. H3, we choose the Lorentzian 
form 



D 2 



D 2 



(5) 



We stress that the hybridization as well as the Coulomb 
interactions to the leads are allowed to be asymmetric, 
which corresponds to a generic setting. Furthermore, we 
do not restrict ourselves to the particle-hole symmetric 
point given by eo = 0. 

We define a£P(tj) = — f= J2k <*( w ~ e fe + Ma) a ia> and 

introduce the vertices 
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for i] = 
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(7) 



The different contributions to the Hamiltonian can then 
be rewritten as 



where we measure the energy uj of the reservoir states 
relative to the chemical potentials p a , 

H D = J2Es\s)(s\, (9) 

s 

with s = 0, 1, and Eq = 0, 2?i = eo respectively, and 

V = rj j duai(u))gi 

+ i J dio J dio'-q-q':a v {Lu')ai{Lu):gi V , (10) 

with the multi-index 1 = r/a containing 77 = ± for cre- 
ation and annihilation operators and the lead index a. 
Repeated incides are assumed to be summed over. We 
consider the case of zero temperature throughout the 
manuscript since temperature is a rather trivial cutoff 
parameter for the RG flow, which at will can be easily 
incorporated in the employed RG formalism. 

III. RG ANALYSIS 

We will study the non-equilibrium properties of the 
IRLM using the real-time renormalization group method 
in frequency spaced (RTRG-FS). The formalism is based 
on a description of the system in Liouville space. The 
density matrix of the full system, p(t), is given by the 
solution of the von Neumann equation 



p{t) = e- iH ^)p{t )e 



iH(t-t ) 



-iL{t-t ) 



(11) 



where L = [H, .] is the Liouvillian acting on usual oper- 
ators in Hilbert space via the commutator. Initially, we 
assume that the density matrix is a product of an arbi- 
trary dot part pz}(io) and grandcanonical distributions 
for the reservoirs, 



p(to) = PD{to) ]Jpre 



(12) 



The object of main interest is the reduced density matrix 
of the dot, which is obtained by tracing out the reservoir 
degrees of freedom 

PD {t) = Tr res p(t) = Tr res e~^^) p D (t ) ]J p a res , 

a 

(13) 

and its Laplace transform (Im 2 > 0) 

p D (z) = / dt e lz{t - 1 ^ p D (t) = Tr res p(t ) . 

J t z-L 

. (14) 

Here the Liouvillian admits the same decomposition as 
(UJ), i.e. L = L res + L*§> + L v with L res = [H res ,.], 



' f du{uj + p a )ai{u)a a (u) , (8) L p ] = \ H D, •], and L v = [V, .]. Using the RTRG-FS we 
J will derive the effective Liouvillian of the quantum dot 
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L e Ji (z) from which the reduced density matrix can be 
calculated via 



Pd{z) = 



PD(k 



(15) 



The stationary reduced density matrix is obtained as 



Pn = lim p D (t) = lim 



TT^Poito). (16) 



The existence of a stationary density matrix was proven 
in Ref. using the RTRG-FS as well as for the Kondo 
model in Ref. [2J using non-equilibrium perturbation the- 
ory to all orders. The matrix elements of the effective Li- 
ouvillian involve the rates for the processes between the 
two eigenstates of the dot, leading to poles of the resol- 
vent (|T5l) at Zp = — iTx and = ±e — iL 2 , where Ti 
corresponds to the charge relaxation rate, T 2 describes 
half of the broadening of the local level, and e is the 
renormalized level position. 

The calculation of the current follows along the same 
lines. The operator for the particle current flowing from 
reservoir 7 to the dot is defined as I 1 = —dN 1 /dt = 
— i[H, iV 7 ], where iV 7 denotes the corresponding parti- 
cle number operator in lead 7. The current in lead 7 
then reads (/ 7 )(i) = Tr£>Tr re;j I 1 p(t). Tracing out the 
reservoir degrees of freedom it can be written as 



(P)(z) = ~iTr D T, 1 (z)p D (z) 



(17) 



in Laplace space, where S 7 (z) denotes the current kernel 
to be derived below. The stationary current is given by 
(7 7 ) s * = -iTrflS 7 (iO+)pg. 

Through L e ^ f {z), the RTRG-FS method provides di- 
rect access to the microscopic cutoff scales. By system- 
atically integrating out the energy scales of the reservoirs 
step by step, a formally exact RG equation can be de- 
rived for L e pf(z) as a function of a flow parameter A, 
where all reservoir energy scales beyond A are included. 
This RG equation is coupled to other RG equations for 
the couplings. Similar schemes can be developed for 
the calculation of the transport current and correlation 
function a 20 ' 22 . All RG equations involve resolvents simi- 
lar to the one occurring in (|15l) , where z is shifted by the 
physical energy scales like the reservoirs electrochemical 
potentials. The cutoff scale is given by the distance to 
resonances, being replaced by the corresponding rate at 
resonance. The microscopic inclusion of decay rates as 
cutoff scales into non-equilibrium RG methods was also 
achieved within flow equation methods 2 ^. 



A. Parametrization and initial conditions 

In Liouville space, defined by the basis (00 11 10 01), 
the bare Liouvillian is given by = [Ho,-], and the 



bare vertices are 



g p(o) = aP 



9v 



for p 
9i for p 



(18) 



and 



r pp'(o) it J for p = 

lr U ; - Opp' \ _ . gxy for p = 



where a + = I and 



10 
10 
0-10 
-1 



(19) 



(20) 



The bare current vertex reads 



p(0) _ 



\r\8 ai pG\ 



p(0) 



For the vertices the following notations are introduced: 



(0) 



Mo) 



£ G P(0) 



^2 G P ir 0) 
p 



together with I^ {0) = E P ( /7 )f ■ We note that ((2TJ) and 
(f2"2")l are related to the commutators and anticommuta- 
tors of ([6]) and ([7]), respectively. In matrix notation, the 
bare Liouvillian and the bare vertices are then given by 

/ 

10 
\ -1 



G (o) 



Mo) 



\P(0) 



p 



(21) 
(22) 



L 



(0) 
D 



G 



(0) 

+0 



f(0) 



1 

-1 








Mo) _ .(o) 



MO) 



G 



MO) 



(0) 



/0 



Vi 1 

/o -1 
1 
110 

Vo 

10 
10 

-110 







1 -1 





1 
1 





/o 

1 

Vo 



with G 



G 



(0) 
-a,+a 



(0) 



-10 
10 





-G 



(0) 



and G 



(0) 



-G 



(0) 



the current vertex we obtain 



Tr£) / 
Trz, I_: 



7(0) 



6 a yt£>(0 1 0) 



'7(0) 



f(0) 

8 aj t^(0 1 



(23) 



For 



(24) 
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Within the RG treatment, the Liouvillian Lu{z) and 
the vertices G\{z;u}\), Gh'(z]uji,uj[) and IJ(z;lui) are 
effective quantities, which obtain an additionai depen- 
dence on the Laplace variable z and depend on frequency 
variables u>i and uj[ (the vertices G x and G^, are only 
needed for the initial setup of the RG flow). In addi- 
tion the current kernel S 7 (z) is generated. As shown 
in App. [A] the dependence of the vertices on the fre- 
quencies lui and can be treated in leading order by 
expanding around uji — lj^ — 0. Therefore, we omit it 
in the following and, furthermore, replace z by its real 
part E = Re{z}. The full z-dependence can be recov- 
ered finally by analytic continuation, which will be done 
in Sec. IIV Bl where we study the time evolution. 

Following Ref. [20l the parametrization of the renor- 
malized quantities follows from charge conservation and 
the following symmetry properties 

Tt d L d (E) = TtdG^E) = Tt d G 1v (E) = , 
L D (E) C = -L D (-E) , £ 7 (£) c - -£ 7 (-£) , 
G^Ef = -a-G- x {-E) , G lv {Ef = G lv {-E) , 

where {A c ) SlSi , , S2S2 , = A*^^^,^ and 1 = -ija. 

As a consequence, the renormalized Liouvillian can be 
written as 



L D {E) 



f -iT + (E) iT-(E) 

iT+(E) -iT-(E) 

e(E) 

V -e{-EY 



(25) 

with T±(E) = T±(~E)*. The renormalized vertices are 
given by 



G+ a {E) — 



G~ a (E) 



and 



/ t a (E) 0\ 

-t a (E) 

0' 

\tl(E) tl(E) 0/ 

/ -t a (-E)* 

t a (-E)* 

tl(-E)* tl(-E)* 

V o ooo 



(26) 



G 



(E) = 



/0 

U a {E) 

\0 







-U a {-E)* 



(27) 



with G^ +a (E) = -G +a ^ a (E). This form of G lv {E) 
holds only in leading order, as higher-order RG contribu- 
tions generate non-zero elements in the upper left 2x2 
block (see Sec. lIII Bl) . while the form ([2"5]) and (f2l)f are re- 
tained to all orders. For the renormalized current vertex 
we obtain the parametrization 



Tr D P +a (E) = -tl(E) (0 10) 



l +a 

n 



as well as for the corresponding current kernel generated 
by the RG flow 



Tr D ^(E)=i(T 1 JE) T 2 (E) 0) , 



(29) 



with T\{E) =m-E)*. 

The bare values, which serve as initial conditions for 
the RG equations, read e(E) — e , T±(E) = 0, t a (E) = 



t 2 a (E) = t 3 a (E) = tg\ U a {E) = U [ a", T\{E) = 0, and 
tl(E) = <5 Q7 ^ 0) . 



(0) 



B. Flow equations 

In this section we summarize the RG equations for the 
renormalized quantities as introduced in the previous sec- 
tion, a detailed derivation is given in App. [A] 

The diagrams taken into account are shown in Fig. 1111 
We consider contributions to the flow of Ld, G\ and G\y 
to lowest order in L ~ t 2 to describe the scaling limit 
and to leading and next-to-leading order in U a to obtain 
exponents up to order 0(U 2 ). Terms of order ~ TU a 
for Gn' are neglected. These would generate nonzero 
elements in the upper left 2x2 block of (|2"T1) . For the 
Liouvillian and the vertices the full independence crucial 
for the time evolution is taken into account. 

Basing on the parametrization of the Liouvillian, the 
current kernel, and the vertices, we introduce the follow- 
ing definitions 



Z(E) .(I - ^e(E) 



T a {E)=2irZ(E + n a )t a (E) 2 , 

r a (E) = Ti(E)-rl(E), 
r' a (E) = ~(T 1 a (E) + T 2 a (E)) , 
r(£) = £r Q (£), r\E) = J2K(E), 

a a 

X (E) = Z(E)(E-e(E)) , 



x'(E)= X (E)~2i l0 Aln 



2A - i x (E) 
A - i x (E) 



(30) 



Tr D Il a (E) = tZ(-E)* (0 1) 



(28) 



where 70 = ^ n ((ji <1 ') 2 , and A is a high-energy cutoff 
which cuts off the Matsubara frequencies of the Fermi 
functions of the reservoirs. Under the RG the cutoff pa- 
rameter A flows from the initial value Ao to zero. The ini- 
tial cutoff is related to the physical reservoir band width 
D by (|A55[) . see App. [A] As shown in App. [A] the flow 
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equations for the effective model parameters read 



C. Analytical solution 



2(UL 0) - 70) 



To 



A-ix'(E + » a ) A-i X '(E + fi a )/2 



(31) 



(1A 



U, 



(0) 



A-i X '{E + n a ) 



IK™ 



2nA-i X '(E + fi a ) 



? a (E) + {E -> -E)* 



f a (E) + (E^-E)* , 



Ta^ = 



(32) 



(33) 



(34) 



/7 (0) 



A + r(^-/i a )-i(^-/x a ) 



The remaining parameters of the Liouvillian and the ver- 
tices are given by 

t ± (e) = ~r(£) ±r'(^) = ±^i /2 (E) , (35) 



^ /3 (£)=i a (£)(l±i7rE/(°)) 
tl{E) = dcnt^E) , 
Z(£)(7 Q (£;) = t/( >. 



(36) 
(37) 
(38) 



As a consequence, it turns out that Z(E)U a (E) is un- 
renormalized up to the second order in the interaction, 
in agreement with previous results^. 

The initial conditions for the RG equations are 



T a (E) = T a (E) 



r (0) 



2 7 r(ii° ) ) 2 , T' a (E) 



and 



X'{E) = E - e + |r(°), where = J2 a ^ ■ For thc 
numerical solution of (|31|) - (f5^t> a discretization in E is 
required, the involved numerical effort is however limited 
due to the fast convergence. 

The RG equations (|3T j) - (|34| reduce to poor man scal- 
ing equations for large A, where all resolvents can be 
replaced by 1/A. In this case similar power laws are ob- 
tained for the stationary current as in Ref. [l6l provided 
that the cutoff parameter is intuitively inserted by hand. 
In contrast, the RG equations derived in this paper reveal 
microscopically the various cutoff parameters. As can be 
seen from (|3T|) - lp3|) . all rates are cut off by the distance to 
resonances, given by x'{E + fi a ). On the other hand, we 
see from (|3"4"|) that the renormalization of the level broad- 
ening, which is contained in the imaginary part of x'(E), 
is cut off by \E — fx a — iT(E — fj, a )\. The RG equations 
presented here go beyond all previous RG analysis for 
the IRLM. Whereas Ref. [H provided a consistent poor 
man scaling analysis without a microscopic derivation of 
the cutoff scales, Refs. [13 and [H showed results from a 
full microscopic non-equilibrium RG analysis, but only 
in leading order in Ua^ for the exponent. 



Within the RTRG-FS approach, the coupled differen- 
tial equations for the flow of the effective system param- 
eters as a function of the infrared cutoff A can be solved 
analytically. The approximate solutions are confirmed by 
numerically integrating the corresponding full RG equa- 
tions (|3"T |l -(|M ll . 

The poor man scaling version of pip , i.e. where the 
resolvents are replaced by 1/A, gives the power-law solu- 
tion 



f Q -> If' (A /A) s ° 
with the exponent 







(0)n2 



(39) 



(40) 



According to ([3"Tt this power law is cut off by x'(E+fJ. a ) 
Therefore, the leading order solution is given by 



r Q (£)^it°) 



A„ 



A-i X '{E + fi a ) 



(41) 



where the exponent is consitently calculated up to 0(U 2 ). 

bmce, for small ui 0) < 1, thc power laws lead only 
to a weak variation, we can use the poor man scaling 
solution (J3SJ) for f a in the other RG equations (|3"2> (1M1). 
and read off the cutoff scale by the remaining resolvents 
in these equations. This gives the following leading order 
solution 



r a (E)~^(t a (E)+f a (-E)*) , 



(42) 



r' a (£)~ (t a ( E )-f a (-E)*), (43) 

X '(E)~E-eo+ l -r e (E), (44) 
with the renormalized level broadening 

An ^ 



r e (E) = J2 T « 



(0) 



A + F(E — fi a ) — i(E — fi a ) 



We note the properties 

T a (E)* =T a (-E) 
T' a (Er=T' a (-E) 

In the limit A — > 0, we obtain 

An 



r e (E) = J2 ri ° } 

(X 

f a (E) = r^ L 



T(E - ii a ) - i(E - fi a ) 
A 



ila 



'jV^E + pb a ) - i(E + [i a - e ) 



(45) 

(46) 
(47) 



(48) 



(49) 



which, together with (l4"2l gives a self-consistent set of 
equations for the determination of T e (E) and T(E). In 
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principle this set can be solved numerically but we will 
provide further analytic evaluations in Section lIVI 

The reduced density matrix Pd{E) = 
(p {E) pi(E) 0) T of the dot in Laplace 
space can be obtained from (| 15[) and ([25]). with 
Lp*(E) = Ld(E)\\ = q. After a straightforward algebra 
we obtain 



Po/i(E) = —Po/i{to) + 



T(E) Po/1 (t Q )-T T (E) 
E(E + iT(E)) 



(50) 



where Po/i(to) are the initial occupation probabilities for 
the dot and T±(E) = V{E)/2±T'{E), according to ([35]). 

Finally using (IT71) and ([2"9")) , the current in Laplace 
space is computed using the density matrix by 



(i a )(E) = -m D i: a (E)p D (E) 

= r 1 a (E)po(E)+rl(E)p 1 (E), 



(51) 



where rl /2 {E) = T' a (E) ± \T a {E), according to 

The stationary probabilities p^^ and the stationary 
current follow frompg^ = lirnE_>.o(— i)Epo/i(E) and 

i st = r^g* + rip? , with r^ = r a (E = o). using 

this gives 



Po/i 


l 

= 2 T 


r 

T : 


j si 


= r' a - 


F' 









where all rates are evaluated at E = 0. 
obtain conservation of probability pg* - 
as current conservation J2 a ^a = 0- 



IV. RESULTS 



(52) 
(53) 

As required, we 
pf = 1 as well 



To determine the level broadening r e (eo), we use (|48|) 
and replace r(eo — — > T(0) = T in this equation by 
using the same arguments as above. In leading order in 
Ua" 1 this gives 



An 



\T + i(fi a -e )| 



(57) 



Neglecting the factor i for the cutoff parameter r e (eo) 
in ([55]) . the self-consistent solution of (|57| and ((55]) is 
approximately 



r,(eo) 



Inserting ([56]) into ([52]) and ([53]) . and using 
find for the stationary dot occupation n st 
stationary current 1^ 



(58) 



581). we 



pf and the 



1 1 \ -s 1 a 
= — I — > — arctan 
o 7T /-^ r 



2 7T 



Ma — CO 

r/2 



(59) 



2r a r ) 3 
f 



j aid an t, ,„ ° — arctan ,„ ) , (60) 



r/2 



r/2 



with Go = \ = ^- m our units. As a consequence, we 
find in leading order the same form as in the noninter- 
acting case (where the result is exact) with Lorentzian 
resonances for the differential conductance at fi a = eo- 
However, the rates T a entering these equations are not 
the bare ones but are strongly renormalized by the in- 
teraction. According to (1551) and ([58"]) they have to be 
determined from the self-consistent equation 



A. Steady-state quantities 

The stationary state is obtained for E = from ([52]) 
and ([53"]) . Inserting the solution ([4*2"]) and (|4"3"1) for T a and 
r^ together with the expression ([4"9"]l for T a , we obtain 
T a = Rcf a and T' a = ^ M Imf a , with 



f - r (0) ( - 

A a — 1 a I 1 



An 



jr e (/Li a ) - i(p a - e ) 



(54) 



Since the cutoff r e (/i a ) is only relevant for \fi a — eo| ~ 
O(r) and since F e (E) varies only weakly as function of 
E, we can replace with good accuracy T e (/j, a ) — » r £ (en) 
in the last equation. Furthermore, neglecting terms with 
higher powers in , we find in leading order 

An 



r Q 



MO) 



l2 F e( £ o) - i(fJ>a ~ £o) 



r' Q ~ -T a arctan ^' 

l e (e )/2 



(55) 
(56) 



1 a — A n. 



An 



||r - i(p, a - e )| 



(61) 



with r = ~^2 a T a . This equation will be further analyzed 
in the next section. In particular, this renormalization 
is responsible for a negative differential conductance at 
large voltage. 

For simplicity, we will restrict ourselves in the following 
mainly to the case of two reservoirs a = L,R with /i^ = 
—fiR — V/2. In this case the dot occupation and the 
current I st = If^ = 



-ijg 4 read 



and 



Go 



T L + V/2 - e 
— arctan ; 

r r/2 



r R v/2 + £o 

— arctan ; 

r r/2 



(62) 



2r L r fl 



arctan 



V/2 - ep 

r/2 



arctan 



U/2 + eo 

r/2 



(63) 
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1. The rates Y a 

As outlined above the rates r a are determined by the 
self-consistent equation (joTj) . We define the cutoff scales 



(64) 



r 

= max{|/i Q - e |, -} . 



From (|6ip we see that T a is renormalized by a power law 
cut off by A^ 



r (o) f'^o X * 
A° 



(65) 



To write this equation in terms of invariant energy scales, 
we introduce the strong coupling scale 



Tk — r|y =£o=0 

and write T a in the form 



I =7T I S 



A:; 



with the independent scales 



r; 



K 



r (o) I ^£_ x 9< 



The scaling limit is defined by — > and Ao — > oo, 
such that Tg remains constant. From (|65|) and (|66l) we 
see that Tr- is determined from the self-consistent equa- 
tion 



(66) 



(67) 



(68) 



(69) 



and remains also constant in the scaling limit. For sym- 
metric Coulomb interactions g a — g, we obtain the solu- 
tion 



T K = r<°> 



lr (0) / 



r (o) 

T K = fk TK ' (?0) 



withrw -E a rL 0) . 

In the special case of two reservoirs a = L/R, we use in 
the following instead of Tg the invariant T K = T K + Tg 
and the asymmetry parameter c 2 = T K /T^. We obtain 



where |/x a — eq I 3> Tjf, charge fluctuations are weak and 
the RG flow of r Q is cut off by the scale \fi a — cq\, which 
describes the distance to the resonance. Close to reso- 
nances, where \fi a — eo| ~ Tk, charge fluctuations are 
strong, and T a is cut off by Tk- Nevertheless, r a is 
bounded by the scale T K for arbitrary system parame- 
ters even for V = eo = 0, leading to finite results for 
all cases. Although there is no rigorous argument why 
our theory should be well-controlled in the presence of a 
single energy scale Tk , we show in the next sections that 
in the scaling limit our results for the charge susceptibil- 
ity and the current are in excellent agreement with exact 

numerical methods, provided that Ua ^ -C 1. This indi- 
cates that strong charge fluctuations are covered by our 
theory. 

Close to resonance, where fi a ~ eo, the rate r Q is loga- 
rithmically enhanced, similar to corresponding logarith- 
mic enhancements for 2-level models with spin fluctua- 
tions (Kondo model), see Ref. [22|. Defining an overall 
cutoff scale by A c = max{A^, A^} and expanding in g a , 
we find close to the resonance 



1 + g a In 



A r 



Mq - e o + *r/2| 



(73) 



In comparison to the Kondo model the IRLM is simpler 
in the sense that the leading order charge fluctuation 
processes provide a unique cutoff scale \fi a — e + iT/2\ for 
the rates. In contrast, for the Kondo model, the distance 
to the resonance as well as the Zeeman splitting itself 
serve as cutoff parameters, such that different logarithms 
can occur for the rates, see Ref. [22| for details. 

The appearance of a negative differential conductance 
in the IRLM at large bias voltages (see Figs. [3] and |4]) 
can be understood^ very easily from the form of the 
rates T Q , while from the numerical or field-theoretical 
computation of the I — V characteristics it is difficult to 
extract the physical mechanisms. In the limit V 3> T, |eo| 
the current (f6"3"|) for two reservoirs reduces to I st (V) ~ 

E^a, and to I st ~ (T ( ° ] T ( ° ] / (T ( ^) 2 )T for symmetric 
Coulomb interactions gi — ga. = .<?• Substituting the 
above expression ([65]) for T being cut off by the voltage, 
and using (|70| . we obtain 



I st (V) 



r (0) r (0) 

L R rp 



T (oy 



K 



Tk 
V 



(74) 



■ T K 



rpR _ 
J w — 



1 



1 + c 2 " ' * l + c 2 
and for symmetric Coulomb interactions 



■ T, 



K 



<0) 
L 
(0) 
R 



(71) 



(72) 



Tk is the energy scale which determines the impor- 
tance of charge fluctuations. Away from resonances, 



leading to a negative differential conductance for re- 
pulsive interactions. The power-law behavior (|74p was 
previously obtained using a variety of other meth- 
ods^ ' 14 i 17 ' 18 i 26 . In contrast, Nishino et alji 2 - find a critical 
value of U = 2 above which negative differential conduc- 
tance appears. For attractive interactions, we obtain a 
power-law increase of the current as a function of voltage 
which is consistent with DMRG results in Ref. 14. How- 
ever, we will show in Section [TV A 31 that this result does 
no longer hold for asymmetric Coulomb interactions. 
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2. Charge susceptibility 

The stationary charge susceptibility \ ( or the static 
capacitance) describes the charge response of the dot due 
to a shift of the level position eo and is defined by 



X 



de 



It can be obtained directly from J5 
level position reads 



1 



X 



2lT „ (Ma 



^o) 2 + (£) 2 



(75) 

and for arbitrary 
(76) 



where we have neglected small corrections from the weak 
dependence of T a on e via the power law (|67|) . For the 
special case of two reservoirs with \xl = —[ir — V/2 and 
for eo = 0, this gives 



r 



XUo=o = - 



7T v 2 + r 2 



(77) 



In particular at V — 0, this result can be compared to 
exact numerical results from NRG, which are shown in 
Fig. H We obtain 



X\t =v=o = 



7rr| £o= y =0 



ttT k 



(78) 



which can be used to define the physical scale Tk even 
away from the scaling limit. For symmetric Coulomb 
interactions <?l = gR = g, we can insert Tk from (|T0[) 
and get 



X|e =V=0 



2 / Ao_V 

7rr(°) ^r(°) y 



(79) 



As can be seen from Fig. [3J the exponent agrees sur- 
prisingly well with the exact numerical result from NRG. 
Since eo = V = is the most critical regime where strong 
charge fluctuations are present, this comparison strongly 
supports that our general solution (|76[) for arbitrary volt- 
age and arbitrary level position is a very good analytical 
approximation to the exact result. 

For | -j ± €q\ ^> r, power laws occur as function of V 
or e . From (j?o) , ([57)1 and (fTTj) we obtain 




For the symmetric case gi = 
2 1 

for V ^> | e 1 , and to 
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TV 
V 



Tk_ 



g this leads to 

2+9 



2+<? 



(81) 



(82) 



U 0.1 



f/ 0i = 0.05/Ti : 





FIG. 2. (Color online) Results for the static susceptibility 
X for the symmetric model with Ijf^ = L^ ' = U 1 - - 1 and 



p(o) _ p(0) _ p(o) at £() _ y _ g. U pp er p ane i : comparison 
of RTRG-FS (solid lines) results, the analytic solution of the 
flow equations (dashed lines), and NRG data (symbols); lower 
panel: exponent j3 = —g/(i + .?)• 



,(0) 



for V < |e |. 

We note that for asymmetric Coulomb interactions the 
general result ([TrJ]) does not exhibit a clear power law if 
the system is coupled to more than one reservoir, neither 
as function of V nor of eo nor of . In this case, a lin- 
ear combination of different power laws is involved which 
does not reveal a clear exponent except for if one of the 
energy scales is much larger than the other two. 



3. Current 

In the case of two reservoirs the stationary current I st 
follows from We take eo > and study the off- 

and on-resonance case separately. A comparison of the 
full numerical solution of the flow equations to the ana- 
lytical results obtained from (|67p and is provided in 
Fig. [3J The excellent agreement shows that for this situ- 
ation already the simplified analytical treatment within 
poor man's scaling yields an accurate description. 

The off-resonance case is defined by \V/2 ± eo| ^> T. 
Using (f6"5)) we obtain 



jst 
jst 



Ll R 



, v 

for — > e 
V for — < e 



G \l-(y/2f 



Inserting for the rates from (|67|) and (|71|) . this gives 



T T K - 



( T K V 



( T K 



( T K V 



1 I T K 
c \JV/2+eo 



(83) 
(84) 

(85) 
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FIG. 4. (Color online) Results for the current /(V) for the 
symmetric model with t ( ° ] = t { £ } = 0.001 and U { ° ] = = 
0.1/-7T off-resonance (solid lines) and on-resonance (dashed 
lines); inset: logarithmic derivative. 



FIG. 3. (Color online) Results for the current I(V) for asym- 

(l±7r/) 0.1/tt with 7(7 = 



metric Coulomb interactions U 



(0) 
L/R 



0.75, 0.5, 0.25, from top to bottom and t£° = t^> = 0.001, 
eo = 0; the numerical solution (solid lines) is compared to 
the analytical result (dashed lines); lower panel: logarithmic 
derivative. 



for V/2 > e , and 



I st = G, 



T 2 
1 K 



(1 + c 2 ) 2 e 2 Q - (V/2) 2 

/ T K \ 9L ( T: 



K 



V (86) 



\\V/2-e \J VI^/2 + eoL 
for V/2 < e . 

From these results one can see in what cases a power 
law can be expected. First, for large voltages V ^> e , 
a power law can only be seen for the symmetric model 
9l = 9b. = g, in which case I st ~ V~ 9 , see Fig. |3l This 
is the same result obtained also in earlier studies&i^ of 
the IRLM. However, in all other cases where gi, ^ gn 
a power law cannot be seen on realistic scales since the 
two terms in the denominator of Eq. (|85[) are typically 
of the same order of magnitude. The scale at which a 
definite power law is recovered is given by the condition 
(T K /V)\ 9L ~ 9R \ /2 < 1 being extremely small for \g L - 
gn] 1. In Fig.|3]the above condition is not met and the 
asymptotic behavior is not observed. Only if in addition 
t° 5l 7^ gn the asymmetry in the bare rates is large 
(c < 1 or c > 1), the power-law behavior of I st (V) is 
recovered (with exponents gh or gn, respectively). 

Interestingly, for V <C eoj a power law also occurs in 
the asymmetric case, since T does not appear in the de- 
nominator of (1561). In this case we obtain 



I st = Go 



(1 + c 2 ) 2 VM 



2+st+gR 



(87) 



i.e. a power law with exponent — (2 + gt + gn) always 
appears as function of the level position eo at fixed volt- 
age. 

In the on-resonance case eo = V/2 the current is given 
by 



2T 



l + c 2 2 C (^ + I(£c)^ 



(88) 



It is important to note that if the level is in resonance 
with one of the reservoirs it is not in resonance with the 
other one. Therefore, at resonance the cutoff scales are 
r for one rate and V for the other. In contrast to the 
off-resonance case, no power law appears even for the 
left-right symmetric model, see Fig. |4j A power law is 
recovered only for unrealistically large V , where the sec- 
ond term in the last denominator of ([55]) can be neglected 
leading to I st ~ V~ 9R . For the symmetric model shown 
in Fig. H the condition (T/V) 9 = 0.01 < 1 is fulfilled 
only for V - 10 30 T K . 

A microscopic determination of the cutoff scales is 
therefore essential to determine the correct on-resonance 
scaling behavior as a function of the voltage, which does 
not simply appear as an additional low-energy cutoff. 
The non-equilibrium physics for the generic situation 
eo = ±V/2 and <?z, ^ gn turns out to be more complex 
and can not be inferred from the linear-response behav- 



4- Conductance 



Another transport property of experimental interest is 
the conductance G — which most vividly features 

the mentioned resonance at eo = ±V/2 as the voltage 
becomes large, see Fig. [5J Analytically, the conductance 
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FIG. 5. (Color online) Conductance G(V) = dl/dV for the FIG. 6. (Color online) Results for the conductance G(eo) for 



symmetric model with t 



(0) 



(0) 



0.001, u 



(0) 



u 



(0) 



0.1 /n and different values of the gate voltage eo = 0, 0.1, 1, 
2, 5, 10 Tk from top to bottom. 



follows from differentiating (|63[) . Neglecting small terms 
from the Independence of the rates r a , we obtain 



G = G 



2T L Tj 

1 r 2 



(i) 2 



(^-eo) 2 + (i) 2 (^ + £0 ) 2 + (i)V 

(89) 



i.e. two Lorentzian resonances at V/2 = ±eo- 
For the off-resonance case V, T <C I eo I we obtain 



r~i ^lT R ~ 

— LtQ 2 — ^0 



(1 + c 2 ) 2 \\e 



T, 



K 



2+gL+gn 



(90) 



in agreement with (|87j) . Results from the solution of the 
full flow equations are shown in Fig. [B] 

On the other hand, as a function of V the current is 
given by 



jst r 2t l t r ( r 
1 =Go ^^{ n - 2 v 



(91) 



for V 3> r, I eo I) where we took into account the first cor- 
rection to the expansion of the arctan- function in (1631) . 
Interestingly, the latter leads to an additional regime 
characterized by a power law independently of the asym- 
metry. Whereas the first term of Eq. does not show 
a power law for asymmetric Coulomb interactions (since 
r = Tl + Tfl appears in the denominator), the second 
term does show a power law because T cancels out. Tak- 
ing the derivative with respect to V, the conductance 
reads 



asymmetric Coulomb interactions U 



(o) 

L/R 



(l±7t/)0.1/7T with 



(0) 



(()) 



7(7 = 0, 0.25, 0.5, 0.75 from top to bottom and t L 
0.001 at V = 0; the obtained power law is characterized by 
the exponent ~ 2{g + 1) in linear order, with g = 2U (0 K 



nates and yields a power law as function of the voltage 



G — Gq 



4£l£r 
V 2 



Gq 



4c 2 



(1 + c 2 ) 2 



Tk 
V 



(93) 



Thus, in contrast to the current, the conductance shows 
always a power law either for large voltage or for large 
level position with the same exponent — (2 + + gn) in 
the range T, |e | < V < T/g. 

However, for V ^ T/g the first contribution in (f9Tj) 
dominates leading to a negative differential conductance. 



B. Time evolution 

The time evolution of the reduced density matrix can 
be obtained directly from (fT5|) via inverse Laplace trans- 
form 



Pn(t) 



i 

2^ 



oo+i0 + 

dz — 

-oo+i0+ Z - 



L^{z) 



Pd(0) , (94) 



G — G( 



f d 2T L T R 

\dv r 



,-21 
V 



2^lTr f 2Y \ Markov approximation L e J^(z) 
V 2 I 



where we have set the initial time to to = 0. We recall 
that we assume the initial density matrix of the full sys- 
tem to be of the product form (fT2"]) . This situation can be 
prepared by setting the couplings between the leads and 
the dot to zero for times t < 0. At t — the couplings 
are suddenly switched on and the system evolves under 
the Hamiltonian ((T|), which results in (f94|) for the time 
evolution of the reduced density matrix. In conventional 



0) 



The derivative in the first term yields a factor 



(92) 
9 IV 



from the weak voltage dependence of the rates T L / R , 
which has to be compared with the factor T/V 2 in the 
second term. Thus, for g <C T/V the second term domi- 



glects the z-dependence of the Liouvillian, which yields 
simple exponential decay towards the stationary reduced 
density matrix. In contrast, we keep the z-dependence 
of the Liouvillian including its branch cuts. The time 
evolution of the current is obtained similarly by inverse 
Laplace transform of (|17[) . This approach has previously 
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been used to study the real-time dynamics of the magne- 
tization and current in the anisotropic Kondo models. 

Specifically, using (|50| we find for the occupation 
n(t) — pi(t) of the dot 

n(t) = (c+cXi) = (1 + J+(t) + J_(i))n(0)- J+(t) , (95) 

where the auxiliary functions J±(t) are defined as 



J± W ^ 2^ 



oo+i0+ 



-00-HO+ 



G?Z 



Z Z + iF(z 



■r±(z) 



(96) 



Here T±(z) = ±T(z) ± T'(z) with T(z) = J2 a r a<» and 
r'( z ) = J2 a ^'a( z ) are obtained from the analytic con- 
tinuation of (g2]), (g31) and (gSl) (see below). Similarly, 
using (|51[) we find for the current in lead 7 



AW 



(97) 



oo+i0 + 



2 

27 



dz e 



-oo+i0+ 
oo+i0+ 



z r(z) 



(ri(z)r_(z) + r2(z)r + (z)) 



"27 



e -«* ri(*)-i*(z) 



dz l, ; 7 .„, 7 . ; (r + (z)-n(o)r(z)) 



-oo+i0+ 



r(z) z + ^r(z) 



where ry 2 (z) = r;(z) ± ir 7 (z). 

In order to evaluate the dot occupation and the current 
in the interacting case we start with the analytic contin- 
uation of Eqs. g2]), d43j), gH) and |49]). We consider the 
case of two reservoirs with — — M-R = V/2 and restrict 
ourselves to the situation of symmetric couplings to the 
leads, i.e. C/f = [#> = and rf = rg> = r(°>/2. 
In this case we can make use of the helpful identity 
T e (-E)* = T e (E), which follows from (05]) and g5]). As 
a result the analytic continuation reads 



r Q (z) 



r(o) 
r(°) 

+ 



E 



An 



T(z - fjt a ) ~ i(z - fi a ) 

Ao 

ir e (z + fi a ) - i(z + fi a - e ) 

Ao 



(98) 



ir e (z 

87rL/(°) 



Ma; 



i(z 



Ma + £0; 

Ao 



(99) 



ir £ (z + ^q,) - i(z + ii a - e ) 

A \ 9 

Aq 



ir e (z- 



Ma 



)-i(z-(JL a + e ) 



(100) 



with ^ = 2i7(°)(l - U^). To proceed we first calculate 
the non-vanishing poles z± = —iT and z± = ±e — ^T e of 
the resolvent l/(z—Lp(z)). Using the parameterization 
([25jl we find — iT(zi) = Z\, e(z+) = z+, and — e(— 2*) = 
z* , which results in 



r(-ir) 



r e (±e- 



r e ±2i(e-e ) . (101) 



Second, we approximate the functions T(z) and T £ (z) in 
the denominators of (f9"5| - (|100j) by their fixed points, i.e. 
we replace 



T(z)^f, r e (z)->r e ±2i(e-e ) 



(102) 



We use the upper (lower) approximation for T e (z) in the 
terms with a singularity at z eo (z w — eo). Inserting 
p02]) we obtain 



r e (*) 
r Q (z) 



^FE 



A" 



TV 
4 



T-i(z- fl a ) 
if e - i(z + /i Q - e) 



(103) 



T, 



if £ - i(z - Ha + e) 



(104) 



r^) = - 



87rJ7(°) 



|f 6 - i(z + n a - e) 

\ 91 



y if £ - i(z - |« Q + e) 
-^- r a( 2 )lnf- -, (105) 



2 1 e 



z(z + /x a - e) 



where have expanded the prefactor in T' a (z) up to 0(1) 
while keeping the exponents unchanged. We note that 
the consistent calculation of 0(U^) corrections of the 
prefactor in T' a (z) requires the analysis of three-loop di- 
agrams which is beyond the scope of our work. Further- 
more we have taken the scaling limit — > 0, Ao — >• 00 
with T]^ 9 = r(°)Ag kept constant. The decay rates f 
and r e and the renormalizcd level position e can be eval- 
uated numerically from Eq. (|101 1) . The physical inter- 
pretation of these quantities is as follows: Whereas T 
describes the charge relaxation processes on the dot and 
thus the relaxation of the diagonal elements of the re- 
duced density matrix with respect to the charge states, 
r e is the broadening of the local level induced by the 
coupling to the leads, i.e. it characterizes the relaxation 
of the off-diagonal elements. Furthermore, the coupling 
to the leads yields a renormalization of the level position 
from the bare value eo to e. For weak Coulomb inter- 
actions, < 0.1, this renormalization is found to be 
small, I e/eo | < 0.01. The approximate analytical expres- 
sions (|103p - (|105p show excellent agreement with the full 
numerical solution of the RG equations (13"TT) - (|3"4")) . 

Inspecting the integral representations and (|9"T)) 
we see that the dominant contributions stem from the 
singularities in the lower half-plane of the involved func- 
tions. We stress that the approximations (|102l) preserve 
this analytic structure, i.e. the poles as well as the po- 
sitions and exponents of the branch cuts from the power 
laws remain unchanged. The integrals (1961) and (|97[) can 
then be treated using standard techniques of contour in- 
tegrations (see App. |B|) . Numerical evaluation yields the 
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FIG. 7. (Color online) Time evolution of the dot occupation 
n(t) for V = eo = 10 Tk and different values of U (0) . The 
initial condition is given by n(0) = 0. We observe oscillating 
behavior at short times, Tut < 5. 



FIG. 8. (Color online) Time evolution of the dot occupation 
n(t) for U m = 0.1/tt, V = 100 Tk and different values of e . 
The initial condition is given by n(0) = 0. 



occupation of the dot n(t) as well as the current I Lit) in 
the left lead shown in Figs. [7] M 

Furthermore, for the long-time behavior off resonance 
(|eo — V/2\ » Tk, 1/t) we are able to derive approximate 
analytical expressions for the dot occupation 27 
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as well as for the current (we assume V ^> Tr-, 1/t in 
addition) 



+^(T jr t)»e- f .*/ a ^ 

Z7T 



r 

'((*4)*) 



> (107) 



with the stationary values given by (j52|) [or ([62]) ] and 
(1531) , respectively. For simplicity we have considered the 
level to be initially empty, n(0) = 0. We note that 
in (|106[) the prefactor of the sine is determined up to 
0(1) while the prefactor of the cosine has been calcu- 
lated consistently up to O(C/ ( - - ) ). The cosine represents 
the leading oscillatory behavior for eo = 0. From Figs. [7]- 
[5] as well as (|106|) and (I107P we observe that the time 
evolution is governed by an exponential decay towards 
the stationary values, characterized by the decay rates 
r and r £ /2. In addition, oscillating terms with explic- 
itly voltage-dependent frequencies e ± V/2 appear, ac- 
companied by an interaction-dependent power-law decay 
~ t 9 ^ 1 . The last result is of particular importance for 
applications in error correction schemes of quantum in- 
formation processing as it violates the standard assump- 
tion^ of a purely exponential decay. The same qualita- 
tive features were observed for the time evolution in the 



anisotropic Kondo model™. We stress that these qualita- 
tive features are independent of the approximations lead- 
ing to (|103p - (|105p as they are completely determined by 
the analytic structure. The imaginary parts of poles and 
branch points lead to exponential decay, their real parts 
yield oscillating behavior, and the integrations along the 
branch cuts result in power laws. 

In Fig. [5] we observe that the current starts at a non- 
zero value. This is due to a non- vanishing displacement 
current^ dn(t)/dt, i.e. the fluctuating number of parti- 
cles on the dot. Specifically, the particle number conser- 
vation in the full system implies 



hit) + I Hit) 



dn(t) 
dt 



(108) 



where lL/B.it) = —dN L /R/dt is the current flowing out 
of lead L/R. The initial condition n(Q) chosen in Fig. [5] 
causes particles to flow from the leads to the dot as soon 
as the couplings t^J R are switched on. These initial cur- 
rents establish on the time scale t ~ 1/D with D de- 
noting the band width, i.e. they start instantaneously 
in the scaling limit. The strong charge fluctuations on 
the dot further result in situations where particles flow 
off the dot into the leads even against the applied bias 
voltage, as can be seen by the appearance of /z,(i) < 
in Fig. |9l A similar displacement current has been ob- 
served by Schmidt et al^S in the transient dynamics of 
the Anderson impurity model, where also the effects of 
different reservoir cut-offs have been investigated. 

In the non-interacting case = the rates and the 
level position are simply given by T — T t — Tk = T*- ', 
and e = eo- The contour integrals l|96p and |97|) can be 
evaluated explicitly (see App. |B|) . resulting in 



n(t) = e~ TKt n{0) + 1 6 ^ + F (t) - F x (t) (109) 



and 



mt)^T K e- TKt 



1 - 2n(0) 



+ ^Y[FoMt)~FoAt)] + ^fF 1 {t), (110) 
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where 7 = -7 and F 0/1 (t) = ^ Q ^b/i,a(*)- Tne nna l 
expression for fo/i,aW is given by the formula 



-T K t/2±T K t/2 
2^ 



± arctan 



Met — £q 
TV/2 



-ImEi ( T-g- + *( £ o - Ma)t 



(111) 



In the stationary limit i — > 00 we recover from (|109[) and 
(jllOp the stationary values ([5^)1 and (fr?5)) . respectively. In 
the opposite limit, at t — + , we use the property (|B9|) 
and observe a non-zero initial value of the current 



Ut = 0+) =T K 



1 - 2ra(0) 



(112) 



i.e. the displacement current discussed above. We note 
that (|109p - (|lll[) contain exponentially decaying terms 
with rates Tk and Trt/2, oscillations with frequencies 
eo i V/2, and power-law behavior ~ l/t. 

Finally, we would like to compare our results for the 
non-interacting model (|109l) - (|lll[) with the literature. 
As is well-known 3 - there exists a mapping (in a certain 
parameter regime) between the resonant level model (and 
thus the anisotropic Kondo model) and the spin-boson 
model (or double-well problem) of dissipative quantum 
mechanics. By studying the time evolution in the latter, 
Lesage and Saleur— showed that generically one has to 
expect relaxation with various decay rates as well as os- 
cillating terms in qualitative agreement with our results. 
Later Anders and Schiller 33 addressed the time evolution 
in the resonant level model. They derived analytic re- 
sults for the dot occupation n(t) in the single-lead model, 
which are identical to the V — > limit of (jl09|) . They 
further considered quite general initial density matrices 
beyond the product form (fT2|) and observed a decay with 
two different relaxation rates accompanied by algebraic 
decay as well as oscillations with frequency eo m this 
more general setting as well. Komnik 3 ^ extended their 
results to the two-lead model and studied the current 
through the system. The results he obtained are simi- 
lar to (|109p and (|110[) . Furthermore, time evolution and 
quench dynamics in the resonant level model have been 
studied in the context of the anisotropic Kondo model at 
the Toulouse poinb 3 -^. 



V. CONCLUSION 

We presented a non-equilibrium RG scheme to study 
transport properties of quantum dots in the regime 
of strong charge fluctuations. We developed a gauge- 
invariant approximation scheme to solve the RG equa- 
tions analytically by expanding all quantities around zero 
Matsubara frequency. We illustrated the approach by 
a minimal and nontrivial model: the IRLM. Whereas 
many previous works treated the problem numerically or 
at the self-dual point, our analytical treatment in the 




FIG. 9. (Color online) Time evolution of the current I hit) in 
the left lead for f/ (0) = 0.1/tt, e = 10 Tk and different values 
of V. The initial condition is given by n(0) = 0. The non-zero 
current at t = is determined by the displacement current 
dn(t)/dt. The current oscillates with frequency ~ eo — V/2; 
note the absence of oscillations on resonance. 



scaling limit and for moderate Coulomb interactions re- 
veals the renormalized tunneling rates r a parametrizing 
the occupation and the current in the same form as the 
nonintcracting solution. The tunneling rates are given 
by a power law cut off by the distance to resonances. At 
resonance the total tunneling rate T — J2 a itself is 
the cutoff scale leading to a self-consistent equation for 
T. We calculated the power-law exponents up to second 
order in the Coulomb interaction and found a very accu- 
rate agreement with NRG. Each r Q has its own power- 
law exponent g a — 2Ua ^ — IlgC^g ) 2 determined by the 

Coulomb interaction between the dot and reservoir 
a. As already pointed out in Ref. [18| it turned out that 
the current does not reveal power laws as function of the 
voltage in the generic case of asymmetric Coulomb inter- 
actions. The reason is that an asymmetry factor T^Tr/V 
occurs, which contains a linear combination of all rates 
in the denominator. In contrast, away from the particle- 
hole symmetric point, we found that power laws occur as 
function of the level position in the generic case, since, for 
large level position, only a factor T^Tr appears for the 
current or the linear conductance. However, the charge 
susceptibility shows a power law neither as function of the 
voltage nor of the level position in the presence of more 
than one reservoir, since it contains a sum of terms, each 
being proportional to the rate r Q . 

Whereas the RTRG-FS scheme is limited to the scal- 
ing limit, the functional RG allows to access the steady 
state for arbitrary tunneling parameters beyond the scal- 
ing limit. As shown in Ref.ll8l both methods provide ex- 
cellent agreement in the scaling limit. The combined use 
of both RG approaches provides hence a complete pic- 
ture of the non-equilibrium physics under consideration. 
As shown in this paper, an advantage of the RTRG-FS 
method is the analytic treatment of Coulomb interactions 
up to next-to-leading order, providing an excellent agree- 
ment of power-law exponents with NRG results. Further- 
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more, the time evolution can be studied with RTRG-FS, 
where we found complex relaxation dynamics similarly to 
previous studies of the dynamics of the non-equilibrium 
Kondo model22. 

The understanding of basic models of spin and charge 
fluctuations opens the way for applications to more com- 
plex quantum dot models. A fundamental issue for the 
future concerns the universality of the effects of strong 
charge fluctuations at resonances found for the IRLM 
where they induce a level broadening and a renormal- 
ization of the tunneling couplings. In particular, in the 
presence of both spin and charge fluctuations, as e.g. in 
the non-equilibrium Anderson model, the level position 
itself becomes renormalized and it is still an open ques- 
tion what the precise line shape of resonances looks like. 
Whereas strong charge fluctuations at resonances seem 
to be described by the RTRG-FS method, an open ques- 
tion remains whether strong spin or orbital fluctuations 
can be covered as well. In both cases, a single energy 
scale dominates the physics and cuts off the RG flow. 
Surprisingly, although no rigorous argument allows the 
truncation of the RG equations in this case, for strong 
charge fluctuations within the IRLM we have shown here 
that our results agree very accurately with exact numer- 
ical methods. Whether such an agreement holds also for 
strong spin and orbital fluctuations will be studied in 
future works. 
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Appendix A: Derivation of the RG equations 

In this Section we report a detailed derivation of the 
flow equations and their evaluation. As outlined in detail 
in Ref. [20|, the RG consists of two steps: The first one is 
a discrete step where the symmetric part of the reservoir 
Fermi function is integrated out, and the second step a 
continuous RG transformation where the Matsubara fre- 
quencies of the reservoir Fermi function are integrated 
out successively (note that at T — the Matsubara fre- 
quencies are continuous). In addition to Ref. ,20 we will 
present here a systematic treatment of the frequency de- 
pendence of the vertices, which can be quite generically 
used for the treatment of strong charge fluctuations. In 
contrast to Rfs. [20j and [22j, where a systematic weak cou- 
pling expansion around the poor man scaling solution 
has been presented, we propose here a systematic expan- 
sion around the point where all Matsubara frequencies 
are set to zero. The dependence on the Laplace variable 
E is fully taken into account, leading to a gauge-invariant 



theory. Many considerations presented here hold in gen- 
eral and can also be used for other models in the charge 
fluctuation regime. In particular, we will make use of 
various generic cancellations of diagrams which simplify 
the RG analysis considerably. 

We consider a model with single and double vertices, 
described by the interaction (flU)) initially. In Liouville 
space the vertices are defined in (TT51) - (|2"2")) . Following 
Ref. HO, we start with the discrete RG step and integrate 
out the symmetric part of the reservoir Fermi function by 
a perturbative treatment. The respective diagrams are 
shown in Fig. [TU]and define the initial values of the sec- 
ond continuous RG procedure. For the initial Liouvillian 
we obtain 

L D (E)=L^- l ^G^- l ^DGf 1 )Of,\ 

-i£G<°Mi' - 4 0) )G$, (Al) 

where we used the short-hand notation Si...„ = E + p,i + 
■ • • + p-n with fli = r)i[i ai . The reservoir band width D 
is related in a certain way to the initial value Ao of the 
continuous RG flow, see Eq. (|A55[) below. 
The initial vertices are given by 

Gl = _ i^Gf - iffi&g (A2) 

and 

The equations for E 7 and J 7 are determined analogously, 
the first vertex just has to be replaced by 17 . We note 

that J7 2 (0) = for our model where double vertices de- 
scribe a Coulomb interaction between the dot and the 
reservoirs, i.e. these processes can not contribute to the 
current vertex. For other models, where double vertices 
describe spin or orbital fluctuations, has to be in- 

cluded as well. 

Inserting the form (|23[) of the initial matrices into 
above expressions and comparing with the parametriza- 
tions (|25p - (|2"9"|) yield the following initial values for the 
continuous RG flow 

r + = r_ = ir(°), r<?> = 2^f , 

e(E) = e (l-^(U^f) + E^(U^f 
t - < (0) f - 5 i (0) 

pl _ i r (0) r 2 _lp(0) 

7 2 7 ' 7 2 7 ' 
U a = . (A4) 



16 



L D 



Gi 



Gn 



1 r 



r 1 



1 r 



1 2 2 2 2 1 



1 2 



2 r 



r 1 



FIG. 10. Diagrams of the discrete RG step. Double vertices 
are represented by dots lying close to each other. The indices 
s/a refer to the symmetric/antisymmetric part of the reservoir 
contraction. 



We proceed with the flow equations for the continuous 
RG procedure. The Laplace variable is decomposed into 
real and imaginary part as z = E + iuj = (E,lu). The 
Liouvillian Ld(E,uj) and the vertices G\{E,uj]lji) and 
Gn' (E, uj; u>i, u>[) acquire an additional dependence on 
the Laplace variable and on Matsubara frequencies uj\ 
and The diagrams taken into account are shown in 
Fig. QTJ We consider contributions to the flow of Lp, 
G\ and G\y to lowest order in T and to next-to-leading 
order in U a to describe the scaling limit and to obtain 
exponents up to order 0(U^). Terms of order ~ TU a 
for Giy are neglected. These would generate nonzero 
elements in the upper left 2x2 block of ([27]) . 
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FIG. 11. Diagrams of the RG equations. Double vertices 
are represented by dots lying close to each other. For G\\i 
diagrams have to be subtracted where the indices 1 and 1' 
are interchanged, provided this gives a new diagram. This 
guarantees the relation Gn< = —Gin. 
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Using the diagrammatic rules developed in Ref. [2(J the RG equations for the Liouvilhan and the vertices read 

- j^L D {E,u) = -i Gi{E, u; A) U(E 1 ,ui + A) G^E^lj + A; -A) 

+ (-if Gi 2 (E, w; A, ui 2 ) U(E 12 ,co + A + cj 2 ) G 2 i(E 12 , to + A + cj 2 ; -w 2 , -A) 
+ (-if G 12 (E) U(E 12 ,lo + A + uj 2 ) G 2 (E 12 ) \J(E x ,lo + A) Gi(E{) 
+ (-if G\(E) U(E u uj + A) G 2 (E\) U(E 12 ,u> + A + ui 2 ) G 21 (E l2 ) 

+ {-tfG 12 (E)U(E 12 ,uj + A + u 2 )G 23 {E 12 )U{E 13 ,u J + A + u 3 )G- 3l (E 13 ) , (A5) 

--^-Gi(£,w;wi) = -iGi 2 (£!,w;wi,A)n(£;i2,w + a;i + A) G 2 (£ 12 , w + a;j + A; -A) 
aA 

-1 G 2 (E, u; A) II(E 2 ,lu + A) G 21 (E 2 ,lo + A; -A, wi) 

+ (-if G 23 (E) U(E 23 ,uj + A + uj 3 ) G 1 (E 23 ) U{E 123 ,uj + ui + A + u 3 ) G 32 (E 123 ) 
+ (-if G 12 (E) U(E 12 ,lu + ui + A) G 3 (E 12 ) n(E 123 ,cu + Wl + A + co 3 ) G 32 (E 123 ) 
+ (-if G 23 (E) U(E 23 ,uj + A + uj 3 ) G 3 (E 23 ) U(E 2 ,u + A) G 2l (E 2 ) 

- (-tf G 23 (E) II(E 23 ,lu + A + uj 3 ) G 31 (E 23 ) IJ(E 12 ,uj + w a + A) G 2 (E l2 ) 

- (-if G 2 (E) U{E 2 ,uj + A) G 13 (E 2 ) U(E 123 , lu + Ul + A + uj 3 ) G 32 (E 123 ) , (A6) 

-^Gn/(.E,w;wiX) = -iGi 2 (E,u- i u) U A) 1 R(Ei 2 ,u) + 0Ji + A) G^ v (E X2 ,oj + u x + A; -A, u[) 

+i G V2 (E, lo; wi, A) n(J5i/ 2 ,cj + cjJ + A) G 5l (£i/ 2 , w + u>[ + A; -A, wi) 

+ (-if 623(E) U(E 23 ,uj + A + uj 3 ) Giv (E 23 ) U(E 1V23 ,uj + ui + wj + A + w 3 ) Sg 3 (£7ii#38) 

- H) 2 G 23 (£) n(£ 23 ,a; + A + w 3 ) G 5l (£; 2 3) U(E 12l lu + ui x + A) G 2V (E 12 ) 
+ (-i f G 23 (E) U(E 23 ,uj + A + lo 3 ) G 3V (E 23 ) U(E V2) u> + u>[ + A) G 21 (E V2 ) 

- (-if G X2 (E) \J(E 12 ,lu + lu x + A) G V3 (E 12 ) \J(E 1V23 , to + w x +cj[ + A + lu 3 ) G 32 (E 1V23 ) 

+ (-if G V2 (E) U(Ey 2 ,uj + lu[ + A) G 13 (E V2 )n(E 1V23 ,u> + lo x + lu[ + A + lj 3 ) G 32 (E 1V23 ) , 

(A7) 



where 



and 



U(E,uj) = 



1 



E 



L d (E,uj) 



(A8) 



G^E) = Gi(£,0;0) , G 1V (E) = G 1V (E,0;0,0) . 

(A9) 

Implicitly, one has to sum over all indices on the r.h.s. 
of (|A5|) - (|A7|) . which do not appear on the l.h.s. In addi- 
tion, one has to perform the integral J Q dui 2 and J Q dui 3 
at all places where the frequencies lu 2 / 3 occur. The RG 
equations for the current kernel and the current vertex 
are analogous to (|A5[> and (|A6[) . respectively, the only 



difference is that the first vertex has to be replaced by 
the current vertex. 

Except for the real part E of the Laplace variable, all 
frequencies are bounded by the cutoff A. Since A — s- 
finally, it is natural to account for the dependence on the 
Matsubara frequencies by expanding the vertices around 
the reference value (IA9|) . where all Matsubara frequen- 
cies are set to zero. Therefore, we have neglected the 
frequency dependence of the vertices in the higher order 
terms of (|A5|) - (| A7|) . The frequency dependence of the 
vertices is calculated in leading order by neglecting the 
frequency dependence of the vertices on the r.h.s. of (|A6[) 
and (|A7I) together with omitting the higher order terms 
in these equations. This gives 



■^{diE^uj-uJt) - diE)} ^ iG 12 (E) (U(E 12 ,u 1 + u 11 +A)-U(E 12 ,A)) G 2 (E 12 ) 

+iG 2 (E) (lJ(E 2 ,ui + A)-U(E 2 ,A)) G 21 (E 2 ) , (A10) 

^{G 1v (E,u 1 ;u} U uj[)-G 1v (E)} c*iG 12 (E) (U(E 12 , u + lo x + A) - U(E 12 , A)) G 2V (E 12 ) 

-iG V2 (E) (II(E 1 , 2 ,uj + lj' 1 +A)-R(E 1 , 2 ,A)) G 21 (E V2 ) . (All) 
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To integrate these equations in leading order we first de- 
fine a function F(E,uj) by 

m(E,ui) = -^F(E,u) . (A12) 
duo 

Neglecting the weak logarithmic A-dependence of 
Ld(E,ui) generated by the RG (note that 11(1?, uj) = 



IIa(£', uj) depends implicitly on A via Ld(E,uj)), we can 
use 

%SI(E, uj + A) ~ -^F{E, uj + A) . (A13) 

Using this in (jAlOj) and (|All|) . and neglecting in addi- 
tion the weak logarithmic A-dependence of the vertices 
generated by RG, we can integrate these equations to 



G 1 (E,uj;uj 1 )~G 1 {E) + G 12 {E) (F(E 12 ,uj + wi + A) - F(E 12 , A)) G 2 (E 12 ) 

+ G 2 (E) (F(E 2 ,uj + A)-F{E 2 ,A)) G 21 (E 2 ) , (A14) 

Gx V {E,u;u u u' x )~G n ,(E) + G 12 (E) (F(E 12 , uj + uj, + A) - F(E 12 , A)) G 2V (E 12 ) 

-G V2 (E) (F(E 1 , 2 ,uj + uj' 1 + A)-F(E V2 ,A))G 21 {E V2 ). (A15) 



To find the RG equations for G\(E), G\y{E) and 



L D {E) = L D {E 7 0) , (A16) 



we set uj — oji = oj[ = in (|A5[) - (|A7|1 . and insert the 
results (|A14|) and (|A15|) for the frequency dependence of 
the vertices in the lowest order terms. Furthermore, the 



I 

frequency integrations can be performed by using 



i [ dujU(E,uj + A) = K{E) , (A17) 
Jo 

with 

K(E) = F(E, 2A) - F(E, A) . (A18) 

Collecting the various terms one finds after some straight- 
forward algebra 



j^L D (E) = iG 1 (E)U(E u A)G 1 (E 1 ) - iG l2 {E)K{E 12 )G 2l (E l2 ) 

- 2i G 12 (E) K(E 12 ) G 23 (E 12 ) K(E 13 ) G m (E 13 ) , (A19) 

-^Gx{E) = iG 12 (E)U(E 12 ,A)G 2 (E 12 ) + i G 2 (E) U(E 2 , A) G 21 {E 2 ) 

+ G 23 (E)n(E 23 ,A + u 3 )G 1 (E 23 )U(E 123 ,A + uj 3 )G i2 (E 123 ) , (A20) 

4tGiv(E) = iG 12 (E)Tl(E 12 ,A)G 21 ,(E 12 ) - i G V2 {E) U(E V2 , A) G 21 (E V2 ) 
aA 

+ G 23 (E)U(E 23 ,A + uj 3 )G 1v {E 23 )U{E 1V23 ,A + uj 3 )G 32 (E 1V23 ) . (A21) 



It turns out that many generic cancellations occur. In 
particular the corrections from the frequency dependence 
of the vertices from the lowest order terms cancel with 
corresponding diagrams in higher orders. For G\(E) 
{Giv(E)) only the first three (two) diagrams of Fig. [TT] 
remain with frequency independent vertices. For Ld(E) 
the cancellation is not complete but the third and fourth 
diagrams cancel against frequency dependent corrections 
of the first two diagrams, whereas the last diagram ob- 
tains a factor 2. These cancellations simplify the RG 



analysis considerably and appear to be a generic model- 
independent feature. 

To calculate the remaining frequency integrations in 
(|A20I) and (|A21|) . and to find explicit representations for 
H(E,uj), F(E) and K(E), we first introduce the spectral 
decomposition of the Liouvillian 

L d {E,uj) = y £,\ i {E,uj)P i {E,w) . ( A22) 

i 

Here, Xi denote the eigenvalues of the Liouvillian and 
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Pi the projectors onto the eigenstates (note that the Li- 
ouvillian is non-hermitian, so that the eigenvalues are 
complex valued, and the right and left eigenvectors are 
not identical). (IA22[) leads to a corresponding spectral 
representation of the resolvent H(E,w), defined in (|A8[) . 
In leading order, we again expand in u) up to first order 
for the eigenvalues Xi(E,w), whereas we neglect the in- 
dependence of the projectors Pi(E,u). This leads to the 
approximation 



Zj(E) 
- iXi(E) 



P(E) 



where we defined the Z-factor 



Zi(E) 



(A23) 



(A24) 



and the distance to the resonant positions 

Xi {E) = Zi(E)(E - Xi(E)) . 

Here, \{E) = \ % (E,0) and P % (E) = P % (E,0). Within 
this approximation we obtain 



(A25) 



F(E,u)~J2 Z iW ~ iXi(E))Pi(E) 



(A26) 



K(E) c 5>(£)ln ( ^.^ffi ^ • ( A27 ) 



The set of RG equations (|A19|I - (|A21|I is thus complete 
and can be solved for a specific model. In particular, it 
turns out that our approximations are gauge-invariant, 
i.e. if all single-particle levels of the dot, all chemical po- 
tentials of the reservoirs and the Laplace variable E are 
shifted by the same amount, all physical observables re- 
main the same. This is only the case if the £ , -dependence 
of all quantities is fully taken into account, an expansion 



around a fixed value of E, like e.g. E — 0, would not 
lead to a gauge-invariant theory. 

We now turn to the evaluation of the RG equations for 
the IRLM. Using the form f|25|) for the Liouvillian, the 
eigenvalues are given by 



X (E) = 0, Xi(E) = —iT(E) , 
X+(E) = e(E), X_(E) = -e(-Ey 

with the corresponding projectors 



(A28) 



Po(E) = 



r_(£) T-{E) 

r+(E) r+(E) o o 



T(E) 



Pi(E) 



P+(E) 



T(E) 






r + (E) 
-r + (E) 
o 










-r_(£0 o o 

T_(E) 





/ 

10 

\o 



P-(E) 



/0 





\ 1 



Using T(E) = r(-E)* 
the ^-functions 



(A29) 

we obtain for the Z-factors and 



= d + ^r(E))- 



Z 1 (E) = Z 1 (-EY 
Xi(E) = -Xi(-E)* = Z^E){E + iT(E)) , 



(A30) 
(A31) 



Z(E) ee Z + (E) = Z_(-E)* = (1 - ^{E))- 1 , (A32) 
X(E) = X+ (E) = X-i-Ef = Z(E)(E - e(E)) . (A33) 

Summing over r\ and taking into account G+ a> - a {E) = 
— G_ Q!+Ct (£'), the flow equations (|A19|) - (|A21|) can be sim- 
plified to 



■^L D (E) = tG +a (E)U(E + » a ,A)G_ a (E + Ma ) + iG_ a (E)Il(E - /j, a ,A)G +a (E - Ma ) 

-2iG +a< _ a (E)K(E)G +a> _ a (E) , (A34) 

-^G +a {E) = iG +at . a (E)U(E,A)G +a (E) - iG +a (E)U(E + fi a ,A)G +at . a (E + fi a ) , (A35) 

(E) — 2G+ Q ' ! - Q ' (E)U(E,A + u> 3 )G (E)U(E,A + uj3)G (E) . (A36) 



The RG equations for the current kernel and vertex read 

^£ 7 (£)=ii; a (£)II(£ + n a + iA)G_ a (E + Mq ) 

+i!l a (E)H(E -fx a + iA)G +a (E - Mq ) , 

A-Il a (E) = -iIl a (E)Tl(E + fx a + iA)G +a< - a {E + » a ) , 

^Il a {E)=iP_ a {E)n{E -fi a + iA)G +a ^ a (E + Mq ) . 

(A37) 



For the derivation of the explicit flow equations for the 
effective parameters as introduced in the matrix repre- 
sentations of the Liouvillian and the vertices we use the 
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following helpful identities: 



GrjaP—rj — PrjGrja — 



G+a,—a'Pl — PlG+ct,—a — , 
P+G- a P- = P-G+ a P+ = , 

G a P i G 3 P j G 4B =0 for i,j = l,±, 
G +a {E)P 1 {E l )G^ = -(t 2 a -t 3 a )(E)t?'(-E'yP_ , 
G_ Q (£)Pi(£i)G +Q , = (t 2 a -t 3 a )(-E)*t<?'(-E')P + , 
G +a ^ a ,(E)P + = P + G +a ^ a (E) = U a (E)P + , 
G +a> _ a ,(E)P_ = P + G +a ^ a (E) = -U a (-E)*P_ , 

(A38) 

and 



G +a (E)P + (Ex)G- a ,{E') = 
t a (E) 



( tli-E'y 
tl(-E'y 






V 

with 

G +ol {E)P + {E x )G^{E' 



tU-E'f 

-tl{-E')* 







(A39) 



-G. a (-EYP-G a (-E'Y • 
(A40) 

Using Eqs. (|A23I) and (|A27|) . the explicit flow equa- 
tions for the rates T±(E) and the level position e(E) are 
determined from the above Eq. (1A34[) for Lo and its 
parametrization (|23|) to 



Z(E + fi a ) 



dA 



A-i X (E + fi a ) 
(E^-EY , 

Zx(E - Ll a ) 



ta(E)ty 3 (-E - jj, a y 



A - iXi(E- Ma) 



(A41) 

t a (E - n a )t' a (E - |U a ) 



with 7(£ l ) 

matrix (|2"B| to 
d 



-2»Z(E) 7 (E)ln 2A ** ( g , (A42) 

£ Q C/a^) 2 and ^(E) = - - 

/i Q J*. Similarly Eq. (|A35[) is evaluated using 



dA 



t a (E) 



Z(E + n a ) 
A-ix{E + ii a ) 



U a {E + fx a )t a (E) , 



which yields 

t a (E) = \[t 2 a {-E - t 3 a (-E - Mq )*] , 

t' a (E)=t 2 a (-E-v a y-t 3 a (-E-n a y = 2mUgh a (E) , 

since the corresponding flow equations have the same 
form and the initial conditions are equal. As a conse- 
quence, introducing the rates T(E) = T + (E) + 



and T'(E) = (T + (E) - T_(E))/2, the flow equations for 
the two rates can be expressed in terms of the single hop- 
ping variable t a (E) . The flow equation for U a is obtained 
by using Eq. (j2"T|) and integrating Eq. (|A36[) over uj 3 



±JJ (F) = 2AZ{E)2 
dA Ua( > (A-i X (E))(2A-i X (E)) 



U a (E) 7 (E) 



(A43) 

Eqs. ([28]) and (f29|) . together with the RG equation for 
the Liouvillian (|A34[) . yield 

^rtl(E) = Z ( E + ^ Ua (E + » a )tZ(E) 
dA A — ix(E + fx a ) 



for the current hopping amplitude, and 
d 

dA 



a-yw -'£ A X^l) g(£)ty ' ( - £ - »° r 

+ {E^-E)* 

for the current rates. Comparing with the equations for 
t a (E) and T±(E) and considering the respective initial 
conditions of Sec. IIIIBl it follows 

tl(E) = S ai t a (E) , r ± (E)=±Y,rl /2 (E) . (A44) 



Summarizing, the flow equations for the effective 
model parameters read 

4rT°W = -2« x Z{E ( + /r ] ^ h a {E) 2 
dA A - i X {E + Ha) 

- (E -> -E)* , 

±r' a (E) = i > Z{E ( + /_r ] J a (E) 2 +(£'—> -EY , 

dA A — tx(E + [i a ) 

± ta {E) = - Z{E ^ U a {E + ^)t a {E) , 
dA A-ix(E + fi a ) 

^re(E) = 2ni^ Z f~/ a) u£h a {E-» a f 
dA A-ix\{E-iia) 



-2iZ(E) j(E) In 



2A - i X (E) 



A-ix(E) 

4rU a (E) = - . 2 ^f^ )2 , U a (E)j(E) , 
dA (A - ix{E))(2A - ix{E)) 

(A45) 

where we introduced the rates T a {E) = T^E) - T 2 a (E) 
and T' a {E) = (1/2) (T^E) + T 2 a {E)), according to the 
definitions (pKl)) . 

We now determine and discuss the equations for the 
Z- factors Z(E) and Z\[E\ and subsequently for x(E) 
and Xi{E). For Z(£) and Zi(£) we find 



We insert the flow equations (|A45|) for £(£*) and T{E) 
and neglect the derivative with respect to E of Zi(E), 
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t a (E), and U a (E) on the right-hand side, as their 
independence is logarithmically weak. This implies 
-^ Xi (E) =s Zi(E)(l - -^K(E)) = 1 and hence 



Z^E 1 - /i Q ) 



(A-iXlf^-ZXa)) 2 

2AZ(£) 



(A-i X (£0)(2A-*x(£0) 



U<?H a (E - 

7(E) 



dA ^ {A - i X (E + Li a )) 2 



(E -» -£)* 



(A47) 



The first terms include an additional factor t 2 a j A and can 
be neglected to leading order, yielding 



± 7(F) „ 2A ZjEfjjE) 

dA 1 ] ~ (A-i X (E))(2A-i X (E)) 



(A48) 



and a constant for Z\{E) ~ 1. The comparison with the 
equation for U a (E) implies 



z(E)u a (E)~uW 



(A49) 



i.e. the product Z (E)U a (E) being unrenormalized. This 
simplifies the flow equation for t a (E) in (|A45|) to 



d , 



u, 



(0) 



A - »x(£7 + A*a) 



f a (.E) 



(A50) 



We can now derive the equation for X (E) = Z(E){E — 
e{E)). The above expressions yield 



with 



d „ Z{E) 



dA 



A - i X i(E - n a ) 



u£h a {E - Mq ) 2 . 

(A54) 

The terms neglected stem from the A-dependence of 
X (E) leading to higher order terms of order "yojx X (E) ~ 
0(U 3 ). The first term of (|A54[) is important since it 
cancels the large term proportional to D in the initial 
condition for e(E), see (|A4|) . To achieve this, the fol- 
lowing relation is needed between the physical reservoir 
band width and the initial cutoff of the RG flow 



A = 



16 In 2 



-D 



(A55) 



Using (|A4[) and neglecting unimportant terms ~ 0(70), 
the initial condition for X '(E) then reads 



X '(E)\ Ao =E-e + i 



r(°) 



(A56) 



The RG equation (|A54p together with the initial con- 
dition (|A56[) lead to the RG equation with the defi- 
nition T a (E) = 2ttZ (E + /i a )t a (E) 2 . Since the first term 
of (|A53[) is of order ~ 70 it can be neglected in all denom- 
inators of (jA45|) and (A - i^E))' 1 ~ (A - ^'(E))- 1 . 
Thus, the first two RG equations of (|A45j) are identical 
to the RG equations (|32l) and (|3"3"|) . Finally, one obtains 
the RG equation (I3T1) for T a (E) if one combines the RG 
equations (|A50[) fo r t a (E ) with the RG equation (|A48[) 
for Z(E) and uses (|A49J) . 



dA 



Z(E)e(E) 



Z(E) 



-270 



A - i X i(E - /jt a ) 
2A - i X (E) 



•In 



Ui°k a (E - Mq ) 2 
AZ(E)e(E) 



A-i x (E) (2A-i X (E))(A-i X (E)) 

(A51) 

with x i(-E) = £ + ir(£) and 70 = £ a (t4 0) ) 2 . From Ec l- 
(|A48[) for Z(J5) the equation for X (E) reads 



dA 



X(E) 



Z(25) 



-270 



i In 



A - «Xi(-B - fl a ) 
2A- i X {E) , 



u^h a {E - 

A X (£) 



A-i X (E) (2A-i x (E))(A-i X (E))_ 

(A52) 



In leading order, this equation can approximately be in- 
tegrated by 

X (E) = 2z 7o A In ^— ^ + X'(£) , (A53) 
A - i X (E) 



Appendix B: Contour integrations for the time 
evolution 



The evaluation of the auxiliary functions J±(t) de- 
fined in (|96[) is performed using standard techniques 
for contour integrations. The integrand has poles at 
z = and z — —iT as well as branch cuts starting at 
z = e ± V/2 - if e /2 and z = -e ± V/2 - if e /2, see 

Fig. m 

Let us first consider the non-interacting case. The so- 
lution of RG equations yields a result for the rates T a 
and T' a , which appears to be exact in the scaling limit. 
In particular, for symmetric coupling we obtain T(z) = 

rf + ii 0) = r<°> ^ t k , r ± (z) = \t k ± £ a r^), 

rl /2 = T' ±\T K , and 



'•' ' IT 



- i(« + /u a - e 
i(z - /i Q + e 



(Bl) 



Using these values we obtain (| 109[) and (|110[) from ([93)1 
and respectively, where the functions F / l a (t) are 



22 



I 



+ 



iOJ )OJ 

I I I 



I + 

iaj ico luj 



V 2 



FIG. 12. (Color online) Analytic structure of the integrand 
of J±(t). All singularities appear in the lower half plane. Red 
dots stand for poles while red solid lines represent branch 
cuts. The pole at z — corresponds to the stationary state. 
The blue dashed and green dotted lines are the original and 
deformed integration contours respectively. 



The branch-cut contribution equals 

jgi(eo-/ia)* 



Fo/i c ; a (t) = 



2ttT k 
VTk , 



d{iy)e 



yi 



iy + i^f- T i^r ~ (eo ~ Ma) 



x— — (— 2ttz) - (e - \i a -> -e + fi a ) 



= — Im ■ 



,i(eo— li a )t- 



Im ■ 



2tt 

o i( e 0-Ma)t- 



I) 



dye 



yi 



2tt 

-T K t/2±T K t/2 



-oo|/Tt + *( £ o - Ma) 
dx e~ x 
x ± 2k* - i(eo - /Lt Q )i 



■ImEi^=F^- + i(eo-Ma)*J) 

(B6) 

where we exploit the analytic continuation of the ex- 
ponential integral function (see 8.212.5 of Ref. HH and 
Ref.EF 



originally defined by 



oo+i0+ 



dze 



2irT K 7_oo+ t o+ z + i^Ti 2 ^ 



We note their important property 



(B2) 



dt 



[FoA*) ~ FiA*)] = T K F lta (t) , (B3) 



which guarantees the current conservation 4;n(t) 



Let us now show that the result of integration in (|B2[) 
leads to Eq. (|lllj) . To this end we deform the contour 
of integration from the real axis to the paths embracing 
the poles and the branch cuts shown in the Fig. [12J The 
positions of the nonzero pole as well as of the branch 
points of T' a (z) (see Eq. (jBll) above) in the integrand 
of (|B2[) are given by the bare values of T — T e = Tk 
and e = eo- We obtain the following contributions to 
Fo/i,^) = Fo/iJt) + F^Jt). The pole contribution 
equals 



F£ a (t) 



r«(o) 



Tj 



K 



, Fl a (t) = r -^±e--« 

J-K 



(B4) 



where 



-T a (-tT K ) = — arctan^-. 



(B5) 



Ei{±z) 



x =F z 



dx, (Re z > 0). 



(B7) 



Combining (|B4[) . (|B5I) and (|B6|) we obtain the formula 
(fTTT|) . We also note that 



F$ 1 %(t = 0+)=-±-Im f° 



dy 



= ± — arctan 

2tt T k /2 



which implies the property 

F 0/1 , Q (t = 0+) = 0. 



-oo|/T Y + *( £ o _ Ma) 

£0 ^ Ma 



(B8) 



(B9) 



In the interacting case, the analytic structure remain 
very similar to that of the non-inteacting case. The main 
difference is contained in the type of branching behavior 
which changes from the logarithmic to the power-law one. 
Additionally, positions of the branch point as well as a 
position of the nonzero pole are shifted to interaction- 
dependent values . 

For |e — V/2\ 3> Tk we can treat all poles and branch 
cuts separately. Thus the evaluation of J± (t) boils down 
to 



/ m- 1 T r(0) 




1 



- y y 

2tt ^ ^ 



-i(a+/3)t 



z + iT{ z y 

- f '' 2 dyeyt 



(BIO) 
-if 



a=±V/2 0=±e 

T±{a + /3-rj + iy) 
a + j3 + iy + iT{a + f3 — rj + iy) 



a + f3 + iy 

-(»?-> -v) 
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Using T±(z) — |r(z) ± T'(z) it can be cast in the form 



2 1 r(o) 



(Bll) 



J b (t) ± JJ(i) , 



where the first two terms equal n st . Then a straightfor- 
ward calculation yields 



2 sin(Trg) e 2 * „ 
■M*) = Z nrr.Ai+ g 1^ 2^ 



-i(a+/3)t 



71- (tT K y+9 



ds 



a=±V/2 /3=±e 



S'J 



A (a+p-ir c /2)t-is , i . 

{tT K y+9 h (-a-2icrt)a ~ l ~ (- 



I 

(to 

I 



(to 

I 



+ 

(tO 

I 



I 

(to 



(to 



+ 

(to 



FIG. 13. (Color online) Integration contour for J±(t) on res- 
onance |e - V/2\ < Tk. 



while Jfc(t) has been evaluated numerically. We note that 
the factor / 3/eju st gives a sign. For V", |e-V/2| > Trt, 1/i 
we obtain (|106l) . Close to resonance, |e — V/2\ <C Tr-, the 
two branch cuts starting at z = ±e=F V/2 — iT e /2 are very 
close to the pole at z = —iT, which leads to a numerical 
instability in the calculation of its residue. Therefore, it 
is advantageous to directly evaluate J±(t) as defined in 
(|96| on the contour shown in Fig. [T3l thereby encircling 
the pole at z = — tT. 

The time evolution of the current given in (1971) can be 



s-2i/3t)s ' (-s-2i(a+p)t)n 1 s 
I 

cast in the form 



is cos(tt 5 ) 



sin 2 (;rg) 
Pi 



where 



/ i (*)=J i (*)+Q-»(o)) j 2 (*)+j 3 w, 



j 2 it) 
j" it) 









'-oc+iO+ 




,.oo+i0+ 




'-oo+i0+ 






2tt\ 





-T' L {z)e-^ 



dz 



T L (z)e 



dz Y'{z)T L (z) 

z z + *r(z) 



(B12) 



(B13) 

(B14) 
(B15) 
(B16) 



The evaluation is analogous to the one for J±(t) above. 
In particular, in the long-time limit Vt, \e — V/2\t 3> 1 
off resonance we find 

w=Fi(z = o) + — (T K ty e {i _ v/2)t - 

J 2 (i)=r L (-if)e- ift , 

j3( ' )= " rW rL(0) + rLHf } e " ft 



(B17) 



where for J 2 and J 3 terms in 0(U^) were neglected. 
The first terms in J 1 and J 3 together yield the stationary 
current fix?)) . 
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